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Abstract 

The SOLA-ECLIPSE Code is being developed to 
enable computational prediction of jet induced 
mixing in cryogenic propellant tanks in a 
low-gravity environment. Velocity fields, 
predicted for scale model tanks, are presented 
which compare favorably with the available 
experimental data. A full scale liquid 
hydrogen tank for a typical Orbit Transfer Vehicle 
is analyzed with the conclusion that coupling 
an axial mixing jet with a thermodynamic vent 
system appears to be a viable concept for the 
control of tank pressure. 


Nomenclature 

C y - turbulence model constant = 0.09 

C e i - turbulence model constant = 1.44 

C e 2 - turbulence model constant = 1.92 

C e 3 - turbulence model constant = 0.79 

G$ - generalized production term 
j - switching integer 

k - turbulence kinetic energy 

P - pressure 

S$ - generalized destruction term 

t - time 

u, v - x and y velocity components respectively 
x, y - spatial coordinates, cartesian or 

cylindrical (x in radial direction) 

- generalized diffusion coefficient 
e - turbulence energy dissipation rate 

ut - "effective" viscosity 

p - density 

(?(< - Prandtl -Schmidt number for k, = 1.0 

o e - Prandtl -Schmidt number for e, o e = 1.3 

$ - generalized variable 


Introduction 

Investigations of liquid propellant behavior 
in a low-gravity environment have been managed and 
actively pursued by the NASA Lewis Research Center 
(LeRC) for over twenty years. The current 
emphasis of the program is directed towards the 
development of technology for the design of fluid 
management systems including liquid acquisition, 
thermal control of cryogenic tankage, and fluid 
transfer. The work being reported is focused on 
problems associated with cryogenic propellant tank 
pressure control . 


* Special Lecturer, Mechanical Engineering. 

** Professor, Mechanical Engineering. 
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The motivation for examining the mixing of 
cryogenic propellants is a desire to efficiently 
employ a Thermodynamic Vent System (TVS) to 
control cryogenic propellant tank pressure. The 
problems associated with a simple direct vent 
system are documented in reference 1. The TVS 
operates by withdrawing liquid from the tank, 
cooling it by expansion through a Joul e-Thomson 
valve, circulating it through a heat exchanger, 
and finally venting it overboard. The propellant 
remaining in the tank is cooled by passing it 
through the heat exchanger and reinjecting it into 
the tank. The efficient application of this 
concept requires that the cooled propellant be 
well mixed with the the other liquid in the tank 
in order to keep all of the liquid below the 
boiling point. 

Scale model experiments using typical 
propellant tank configurations have been performed 
in the LeRC Zero-G Facility. Reference 2 reports 
the results of experiments which used ethanol in 
plexiglass tanks to simulate the behavior of the 
propellants. Dyed ethanol was injected to 
simulate a mixing jet and high-speed photography 
was used to record the mixing. Ethanol was chosen 
as the experimental fluid because the use of 
actual cryogenic propellants was impractical. The 
experiments were performed in 10 cm diameter tanks 
because only five seconds of free fall is 
available for the experiment. In this size tank, 
approximately two seconds is required for the 
liquid-vapor interface to reach a zero-g 
configuration, leaving three seconds of 
experimental time to examine the liquid jet mixing 
phenomena. A dyed liquid jet was used due to the 
difficulty of establishing and measuring 
temperature gradients in small scale tanks with 
limited test time. The experimental results were 
used to develop flow characterization and mixing 
parameters which can be used to determine mixer 
requirements for full scale tanks. 

Several factors provide motivation for 
developing the capability of computationally 
predicting these flows. A computational analysis 
is not restricted to small tanks, short times, or 
simulated propellants. It can predict both the 
velocity and temperature fields due to the 
injection of a cool jet. Further, it will provide 
a more economical tool for investigating design 
and operational paramaters. 

The SOLA-ECLIPSE Code (Energy Calculations 
for LTcju id Propellants in a Space Envi ronment ) is 
being developed to examine jet induced mixing of 
cryogenic propellants in a low-g environment. The 
starting point for its development was the NASA 
S0LA-V0F Code^ which was written to analyze the 
draining of cryogenic propellants from tanks in 
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low-g. The following sections present a brief 
description of the development of SOLA-ECLIPSE, 
the results which have been obtained by modelling 
the drop tower experiments 4 , and preliminary 
results for typical Orbit Transfer Vehicle (OTV) 
liquid hydrogen tanks. 


Outline of NASA-SOLA-VOF 

The NASA SOLA-VOF Code is a derivative of the 
SOLA family of codes which use the SOLA algorithm 
to solve the two-dimensional (Cartesian or 
cylindrical) Navier-Stokes equations on an 
Eulerian mesh of rectangular computational cells. 
The cell size may be varied by row and by column. 
The code solves the laminar hydrodynamic problem 
for the primitive variables P, u, and v. A 
staggered grid is used with velocities evaluated 
at cell faces and pressure defined at cell 
centers. 

The VOF (Volume-of -Fluid) concept is based on 
a function, F, which has a value of one at any 
point in space occupied by fluid and a value of 
zero at a point which is not occupied by fluid. 
With F so defined, the fractional volume of a cell 
which contains fluid is the average value of F for 
that cell. Cells through which a free surface 
passes are identified by values of F between zero 
and one. The free surface location and shape is 
determined by evaluating F in the cell being 
considered and in its neighbors. The curvature of 
the surface is computed, and surface tension 
effects are modeled by imposing an appropriate 
pressure at the free surface. 


Turbulence Model 

The NASA SOLA-VOF code solves the laminar 
hydrodynamic problem. Since the goal of injecting 
the jet into the tank is to induce mixing, it is 
anticipated that a turbulent jet will be used due 
to its superior mixing characteristics as compared 
to a laminar jet. The following paragraphs 
describe the turbulence model which was installed 
in SOLA-ECLIPSE. 

The k-e turbulence model was chosen as the 
turbulence model for the code because it is the 
simplest model which has a good probability of 
succcessful ly modeling the recirculating flow 
which is induced by the mixing jet. It is the 
most widely used two equation turbulence model and 
its success with flows of this type is well 
documented^ . This model results in the 
computation of an "effective viscosity" which is 
then used in the momentum equations. The two 
differential equations of the k-e model, and the 
equation to compute the "effective" viscosity, can 
be written in the following form for two 
dimensional incompressible flow: 




where: $ = k, e 


.* _ 0 - Cartesian 


J = 


1 - cylindrical 


M* = C p - 
t u e 


( 2 ) 


k-equations 
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where: 


S k = pe 

(5) 

e-equations 
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Although this set of expressions and values has 
produced good results for plane jets and for wall 
boundary layers, it overpredicts the spread rate 
of axisymmetric jets by approximately forty 
percent. Several functional dependencies for the 
constants have been proposed to cure this 
problem^*?*^ # but all the proposed relationships 
require computations involving jet width and/or 
centerline velocity gradient which would greatly 
complicate computations. Reference 9 proposes 
that the existing constants and expressions be 
retained with the inclusion of the following 
expression as an additional source term in the 
e-equation. 


C 


e3 


k£.i/ 9 V 3u 
e 3x ” 3y 



( 10 ) 


This term accounts for the vortex stretching which 
occurs in axisymmetric jets but is not present in 
true 2-D flows. The value of C e 3 was chosen by 
the author of reference 9 to produce the correct 
spread rate for axisymmetric jets. 

The question of appropriate solid boundary 
conditions for this turbulence model has been 
addressed by many workers. A law-of-the-wal 1 
boundary condition has been installed in 
SOLA-ECLIPSE in order to avoid the cell packing 
required to resolve the steep gradients near a 
wall. This approach replaces the computations of 
the gradients by approriate expressions derived 
from a law-of-the-wal 1 analysis. The actual 
expressions used in the code are developed in 
reference 10. 

The free surface boundary conditions for the 
turbulence model were chosen to be as simple as 
possible since the free surface activity is not 
the focus of this study. Since turbulence aids 
the mixing process, underpredicting turbulence 
will result in conservative results. This can be 
achieved by underpredicting the turbulence kinetic 
energy and overpredicting the turbulence energy 
dissipation rate. This is accomplished in the 
code by setting k equal to zero in a surface cell 
and setting e equal to the value of e in the 
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adjacent full cell. Although no claim is made 
that the code accurately models turbulence effects 
at the free surface, the accuracy of the results 
obtained for the bulk fluid behavior should not be 
significantly affected. 


Computational Results for Model Tanks 

As previously described, experiments have 
been performed to examine jet induced mixing in 
zero-g by injecting a dyed jet into partially 
filled scale model tanks. The first cases to be 
modeled with SOLA-ECLIPSE were replications of 
experimental conditions. The following paragraphs 
discuss the predicted velocity fields for one 
laminar jet and for one turbulent jet for which 
experimental data is available. 

The computational mesh used to analyze the 
scale model tanks is shown in Figure 1 and the 
resulting code approximation to the tank is shown 
in Figure 2, The variable mesh capability of 
SOLA-ECLIPSE has been used to refine the mesh in 
regions which are expected to contain sharp 
gradients. The vector field displays show both 
the true curve of the tank head and the stairstep 
approximation to it. For both of the cases being 
presented, the upper head has been removed since 
no flow occurs in this region. A scaling vector 
is supplied so that vector magnitudes can be 
interpreted. The experimental tanks are 10cm in 
diameter and the scaling velocity is given in 
cm/sec. Velocities smaller than two percent of 
the jet velocity are not displayed. Fluid is 
injected through the 0.4 cm diameter pipe which 
penetrates the lower head and is simultaneously 
withdrawn through the annulus surrounding the 
pipe. The heading includes the elapsed time from 
jet initiation (seconds) and the number of cycles 
(time steps) required by the code to reach this 
point in the transient. The detail of the free 
surface shape very close to the wall is 
misrepresented by the graphics but is accurately 
represented in the calculation procedure. The two 
scale model cases have been run to an elapsed time 
of three seconds which was the maximum time 
available for mixing in the LeRC Zero-6 Facility. 

The development of the laminar jet is 
depicted at one second intervals in Figures 3, 4, 
and 5. The highly curved free surface is a 
reflection of the zero-g environment where surface 
tension is an important force. The jet Reynolds 
number for this flow is 450 which is well below 
the transition value. The experimental data 4 
revealed no measurable spreading of this jet. 
Comparison between the predicted velocity field 
and the experimental results is clouded by the 
inherent differences in information being 
evaluated. The experimental data consists of 
pictures of dye spreading which really represents 
the concentration field being produced. Hence, 
the picture only reveals the current location of 
fluid which was injected. In contrast, the 
predicted velocity field shows both the motion of 
fluid which originated in the pipe and the induced 
flow of the fluid which initially resided in the 
tank. It is therefore necessary to estimate where 
the edge of the jet indicated by dye diffusion 
would lie in the predicted velocity field. 
Examination of the predicted velocity field shows 
very little induced flow and leads to the 
inference of very little jet spreading. Figure 3 
shows a slight mushrooming of the head of the jet. 


This is characteristic of the the jet development 
prior to reaching the free surface and is in 
agreement with the experimental data. Figures 4 
and 5 show a slight bulge in the free surface at 
the centerline of the tank. Although the bulge is 
more gradual and of a smaller height than the 
experimental results, the difference is not large 
and it does not change the fundamental nature of 
the flow. Examination of the velocity field after 
three seconds of jet operation indicates 
approximately the same spread of jet fluid along 
the surface toward the tank walls as was found in 
the experiment. 

The turbulent jet case modeled with the code 
used the same physical configuration as the 
laminar case. The jet velocity is 50.1 cm/sec 
which corresponds to a jet Reynolds number of 
1320. The predicted velocity fields are shown in 
Figures 6, 7, and 8. An acceleration of 
0.0022 g's was specified to match the experimental 
conditions. As can be seen by comparison whith 
the laminar case which was in a zero-g 
environment, even a very small acceleration will 
have a significant impact on the curvature of the 
free surface. Examination of the predicted 
velocity fields reveals several interesting 
features. First, the induced velocity field is 
much more extensive than was predicted for the 
laminar jet. Although the edge of the jet is 
again difficult to determine, a ten degree spread 
rate appears reasonable based on the velocities 
displayed. This is in good agreement with both 
theoretical predictions for steady axisymmetric 
jets and with the experimental results 4 . The 
centerline bulge is again underpredicted but the 
character of the experimental and predicted flows 
is the same. The experiment shows dye spreading 
all the way to the tank walls and starting to flow 
back toward the tank head. The field displayed in 
Figure 8 appears to be in agreement with this 
level of development. Figure 6 shows considerable 
free surface activity, even before the jet reaches 
the free surface. In order to compare the 
predicted level of activity to the experimentally 
observed activity, a movie film was made by 
sequentially displaying predicted velocity fields. 
Comparison of the computer generated film and the 
experimental data film shows comparable levels of 
free surface activity. 


Orbit Transfer Vehicle Tank 

The first full scale tank to be modeled was a 
liquid hydrogen tank for a typical Manned Orbit 
Transfer Vehicle (0TV) 11 . The tank is slightly 
more than two meters in radius and the cylindrical 
section has a length of almost six meters. 

Although the proposed tank has elliptical heads, 
the code does not currently have the capability of 
modeling elliptical heads therefore hemispherical 
heads were used for the computational model. The 
computational mesh used for analysis of this case 
is shown in Figure 9. The tank was assumed to be 
fifty percent full in a zero-g environment which 
resulted in the initial conditions depicted in 
Figure 10. The correlations developed in 
reference 4 were used to select a jet velocity 
which would result in a deformation of the free 
surface but would not result in a large geyser. A 
tank-to-jet radius ratio of 22:1 was chosen to 
approximate the ratio of the experimantal tanks of 
the previous section. A velocity of 4.51 cm/sec 
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was chosen which corresponds to a jet Reynolds 
number of 50,000* The code was run to an elapsed 
time of two hours. It was felt that this was 
sufficient time to determine the feasability of 
using a mixing jet in conjunction with a TVS to 
cool the bulk of the propellant in this 
configuration. 

The development of the velocity field 
predicted by SOLA-ECLIPSE for this case is shown 
in Figures 11, 12, 13, and 14. After five minutes 
of operation, the jet has reached the free surface 
and is spreading radially. The predicted field 
after twenty minutes of jet operation shows that 
at least a small velocity has been induced in most 
of the fluid. Comparison of the predicted fields 
after one and two hours indicates that the 
velocity field is fully developed after one hour. 
Almost all of the fluid in the tank has been set 
in motion. After one hour, the large 
recirculation region which was anticipated has 
appeared and extends from the free surface down 
into the upper part of the lower head. 

These results can be compared to the 
correlations for geyser height and for mixing time 
developed in reference 2. The correlation 
predicts a geyser height of 0.4 tank radii whereas 
SOLA-ECLIPSE predicts a geyser height of 
approximately 0.1 tank radii. Although agreement 
on geyser height seems poor, the character of the 
mixing is not significantly different for the two 
predictions. In both analyses the jet does not 
pierce the free surface but is deflected toward 
the tank walls. The mixing time correlation 
predicts the time required to completely mix the 
fluid in the tank. For the situation being 
analyzed, the correlation predicts a mixing time 
of 3700 seconds. Since the velocity field 
predicted by SOLA-ECLIPSE appears to have reached 
steady state by this time, the conclusion of 
thorough mixing seems reasonable. 

A word of caution must be voiced concerning 
the conclusions which can be drawn from these 
results. The purpose of the jet is to disperse 
the cooling effect of a TVS in order to subcool 
the bulk of the propellant in the tank. The goal 
is therefore not to simply induce propellant 
motion, but to suppress propellant temperature 
throughout the tank. Although the velocity field 
may appear to be fully developed, this does not 
insure that the propellant temperature is 
significantly reduced at this time in the 
transient. The extent of the induced motion does 
indicate that convection should be a fairly active 
process and that cooling of most of the field 
should be possible. Work is currently in progress 
to add a thermal energy equation to SOLA-ECLIPSE 
and thereby directly predict the relationship 
between fluid temperature and time. 


Conclusions 

The following conclusions have been drawn 
based on the computational analysis of flows in 
scale model tanks. SOLA-ECLIPSE is capable of 
predicting the small spread rate of laminar jets 
without introducing excessive artificial 
diffusion. The code is capable of modeling flows 
which involve some agitation of the free surface. 
The turbulence model in SOLA-ECLIPSE predicts 
appropriate spread rates for axi symmetric 
turbulent jets. Therefore, the code is judged to 
be a viable tool for studying jet induced mixing 


in cryogenic propellant tanks in a low-gravity 
environment. 

The use of an axial jet to disperse the 
cooling effects of a thermodynamic vent system, 
and to thereby control tank pressure, appears to 
be a viable concept. Although actual temperature 
distributions have not been computed, the velocity 
field induced by the jet extends throughout the 
tank indicating that the fluid will become well 
mixed. Clearly, more extensive analyses are 
required before complete proof of concept can be 
claimed. 
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Figure 1 . - Computational mesh used for analysis of scale 
model tanks. 
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Figure 2 . - Computer code approximation of scale model tank in- 
cluding laminar jet initial conditions. 





Figure 3. - Calculated velocity field with laminar jet after 1 sec 
(1992 computational cycles). 


Figure 4 - Calculated velocity field with laminar jet after 2 sec 
(3894 computational cycles). 





Figure 5. - Calculated velocity field with laminar jet after 3 sec 
(5795 computational cycles). 


Figure 6. - Calculated velocity field with turbulent jet after 1 sec 
(2086 computational cycles). 






Figure 8. - Calculated velocity field with turbulent jet after 3 sec 
(5927 computational cycles). 


Figure 7. - Calculated velocity field with turbulent jet after 2 sec 
(4010 computational cycles). 
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Figure 9. - Computational mesh used for analysis of full scale 
tanks. 


Figure 10. - Computer code approximation of full scale liquid 
hydrogen tank including turbulent jet initial conditions. 





Figure 11. - Calculated velocity field for OTV liquid hydrogen 
tank with turbulent jet after 5 min (1315 computational 
cycles). 


Figure 12. - Calculated velocity field for OTV liquid hydrogen 
tank with turbulent jet after 20 min (2952 computational 
cycles). 
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Figure 13. - Calculated velocity field for OTV liquid hydrogen 
tank with turbulent jet after 1 hr (7320 computational cycles). 


Figure 14. - Calculated velocity field for OTV liquid hydrogen 
tank with turbulent jet after 2 hr (14244 computational cycles). 
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